Self-averaging of Random and Thermally disordered diluted Ising systems 
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Self-averaging of singular thermodinamic quatities at criticality for randomly and thermaly di- 
luted three dimensional Ising systems has been studied by the Monte Carlo approach. Substantially 
improved self- averaging is obtained for critically clustered (critically thermally diluted) vacancy 
distributions in comparison with the observed self-averaging for purely random diluted distribu- 
tions. Critically thermal dilution, leading to maximum relative self-averaging, corresponds to the 
case when the characteristic vacancy ordering temperature (9) is made equal to the magnetic crit- 
ical temperature for the pure 3D Ising system (T,? D ).For the case of a high ordering temperature 
(6 >> T^ D ), the self-averaging obtained is comparable to that in a randomly diluted system. 



PACS numbers: 05.50.+q, 75.10.Nr, 75.40.Mg, 75.50.Lk 
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Systems with quenched randomness have been studied 
intensively for several decades jjj. One of the first re- 
sults was the establishment of Harris criterion |2] , which 
predicts that a weak dilution does not change the char- 
acter of the critical behavior near second order phase 
transitions for sustems of dimension d with specific heat 
exponent lower than zero (so called P systems), a pure < 
==> Vpure > 2/d, in the nonrandom case. This crite- 
rion has been supported by several renormalization group 
(RG) analyses Hf|,U, an d by scaling analysis M. It was 
shown to hold also with strong dilution by Chayes et al. 
[j/J. For a pure > (called R systems), for example the 
Ising 3D case, the system fixed point flows from a pure 
(undiluted) fixed point towards a new stable fixed point 
at which a pure < |g|,|f7g]. 

For a random hypercubic sample of linear dimension 
L and a number of sites N = L d , any observable singu- 
lar property X presents different values for the different 
realizations of randomness corresponding to the same di- 
lution. This means that X behaves as a atochastic vari- 
able with average [X], variance (AX) 2 and a normalized 
square width Rx = (AX) 2 /[X] 2 . A system is said to 
exhibit self averaging (SA) if Rx — * as L — > oo. If the 
system is away from criticality, L >> £, being £ the cor- 
relation lenght, the central limit theorem indicates that 
strong SA must be expected. Hovewer, the behavior of 
a ferromagnet at criticality, where £ >> L, is not so ob- 
vious. This point has been studied recently. Wiseman 
and Domany (WD) investigated the self-averaging of di- 
luted ferromagnets at criticality by means of finite-size 
scaling calculations || , concluding weak S A for both the 
P and R cases. In contrast Aharony and Harris (AH), 
using a renormalization group analysis in d = 4 — e di- 
mensions, proved the expectation of a rigorous absence 
of self-averaging in critically random ferromagnets p0[ . 
More recently, WD used Monte Carlo simulations to 
check the lack of self-averaging in critically disordered 
magnetic systems O] . The absence of self-averaging was 



confirmed. The source of discrepancy with their previ- 
ous scaling analysis was attributed to the particular size 
dependence of the distribution of pseudocritical temper- 
atures used in their work. 

Quenched randomness has been investigated basically 
under two different constrains, a grand-canonical con- 
strain (average density of occupied spin states (p) fixed) 
and a canonical constrain (total number of occupied sites 
(c), constant). The appicability of the result obtined 
by Chayes et al. also to canonical ensembles, imply- 
ing universality between both kinds of constrains, has 
been recently investigated [ Qp^ ]. Monte Carlo simula- 
tions |H],[n| indicate universality for Rx in the grand- 
canonical constrain ensembles with different concentra- 
tions (p). Universality between results obtained in the 
canonical and the grand-canonical constrains at very 
large values of L is implied by the work of Aharony et al. 

In all cases previously investigated frozen disorder was 
always produced in a random way, that is, vacancies 
were distributed throughout the lattice randomly. Real 
systems, however, can be realized in other kinds of va- 
cancy distributions. In particular, dilution in a ferromag- 
netic lattice can be produced by an equilibrium thermal 
order- disorder distribution of vacancies, governed by a 
characteristic ordering temperature (9). If this order- 
ing temperature 6 is high enough, the equilibrium ther- 
mal disorder will be similar to the random disorder of 
previous investigations. However, if 9 happens to coin- 
cide with the characteristic magnetic critical temperature 
(T c = 4.511617) of the undiluted system, new possibili- 
ties are open. In the present work we study the effects on 
the magnetization self-averaging in diluted ferromagnets 
obtained by thermal vacancy distributions in three di- 
mensional Ising systems using the Monte Carlo approach. 
We will study the self-averaging or lack of it in these sys- 
tems and will compare the results with those obtained 
with random vacancy distributions for the same concen- 
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FIG. 1. Scattered plot of the concentration of spins (c) 
versus the normalized magnetization per spin M/[M] ([M] 
is the media of the magnetization in all the realizations: 
[M] random » [M] thermal » 0.15). Results are shown for the 
random case with p = 0.5 and c > 0.5 (open circles) and 
the hipercritical thermal case with 6 = 1000 » T'i D (black 
circles). Dispersion of the magnetization and dispersion of 
the concentration are similar in both cases. Inset shows the 
histograms for the values of the normalized magnetization in 
both cases. The upper case corresponds to random dilution 
and the lower to thermal dilution. Note the equivalence of the 
shapes resulting in similar values for the normalized square 
widths. 

at the critical temperature for a dilution equal to 
p = 0.5, T c {p = 0.5) in order to compare the normal- 
ized square width in both cases. In this case we find ex- 
tremely low values of the normalized square width for the 
thermally diluted set of samples in the " grand-canonical" 
case, R t ^ I ermal « 10~ 3 , which is around 60 times smaller 
than the width found for the random case. (See the width 
differences in the histograms on Fig. 2). 

To study the self-averaging evolution we must consider 
the behavior of the normalized square width Rm with the 
length L of the system. Consequently we have diluted the 
systems thermally at criticality (9 = T^ D ) for different 
values of L. Of course, in all this realizations for critical 
thermal cases, the corresponding value of p depends on 
the value of L, but tends always to 0.5 as L increases 
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FIG. 3. Scattered plot of the concentration of spins (c) ver- 
sus the normalized magnetization per spin M/[M] for critical 
thermal dilution at different lengths (L), L = 20 (black cir- 
cles), L — 30 (open circles) and L = 40 (black squares). We 
find [M] thermal ps 0.9 for L = 20, [M] thermal « 0.7 for L = 30 
and [M] thermal « 0.6 for L = 40. Note the decrease in the 
dispersion of c as L increases. This behavior means that, as L 
increases, the system tends to a constant value of c = 0.5 cor- 
responding to the canonical constrain. Inset shows the same 
scattered plot but now for the random case with p = 0.5 and 
L — 20 (black circles) and L = 40 (open circles). 
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"grand-canonical" samples with L increasing, show sub- 
stantial self-averaging improvement of the magnetization 
in three dimensional highly diluted (p = 0.5) Ising ferro- 
magnets. 
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FIG. 4. Log-Log plot of the normalized squared width of 
the magnetization Rm vs. the lenght L of our system for the 
"grand canonical" case. Results are shown for the random 
dilution with p — 0.5 and c > 0.5 (open circles) and for the 
critically thermally dilution (black circles). 

In conclusion, our Monte Carlo results in critically 
thermally diluted (as opposed to randomly diluted) 
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